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Abstract 

The presence of a particle continuum, both of a resonant and non-resonant character, 
can significantly impact spectroscopic properties of weakly bound nuclei and excited 
nuclear states close to, and above, the particle emission threshold. In the framework 
of the continuum shell model in the complex momentum-plane, the so-called Gamow 
Shell Model, we discuss salient effects of the continuum coupling on the one-neutron 
overlap integrals and the associated spectroscopic factors in neutron-rich helium and 
oxygen nuclei. In particular, we demonstrate a characteristic near-threshold energy 
dependence of the spectroscopic factors for different ^-waves. We show also that the 
realistic radial overlap functions, which are needed for the description of transfer 
reactions, can be generated by single-particle wave functions of the appropriately 
chosen complex potential. 
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Overlap integrals, Exotic Nuclei 
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1 Introduction 



Today, much interest in various fields of physics is devoted to the study of 
small, open quantum systems (OQS), whose properties are profoundly affected 



Preprint submitted to Elsevier 



1 February 2008 



by the environment, i.e., continuum of decay channels. Although every finite 
fermion system has its own characteristic features, resonance phenomena are 
generic; they are great interdisciplinary unifiers. Many of these phenomena 
have been originally studied in nuclear reactions, and now they are explored 
in molecules in strong external fields, quantum dots and wires, and other solid 
state micro- devices, crystals in laser fields, and microwave cavities. 

In the field of nuclear physics, a growing interest in the theory of OQSs is 
associated with experimental efforts in producing weakly bound and unbound 
nuclei close to the particle drip lines, and studying structures and reactions 
with those exotic systems. In this context, the major challenge for nuclear 
theory is to develop theories that would allow to understand properties of 
those exotic physical systems possessing new and different properties [1,2,3]. 
To this end, a unification of structure and reaction aspects of weakly bound 
or unbound nuclei, based on the OQS formalism, is called for. 

The nuclear shell model (SM) is the cornerstone of our understanding of nu- 
clei. In its standard realization [2,4], the SM assumes that the many-nucleon 
system is perfectly isolated from an external environment of scattering states 
and decay channels. The validity of such a closed quantum system (CQS) 
framework is sometimes justified by relatively high one-particle (neutron or 
proton) separation energies in nuclei close to the valley of beta stability. How- 
ever, weakly bound or unbound nuclear states cannot be treated in a CQS 
formalism. A consistent description of the interplay between scattering states, 
resonances, and bound states in the many-body wave function requires an 
OQS formulation (see Ref. [3] and references quoted therein). Properties of 
unbound states lying above the particle (or cluster) threshold directly impact 
the continuum structure. Coupling to the particle continuum is also impor- 
tant for weakly bound states, such as halos. A classic example of a threshold 
effect is the Thomas- Ehrman shift [5], which manifests itself in the striking 
asymmetry in the energy spectra between mirror nuclei. (See Ref. [6] for more 
discussion and examples related to this point.) 

In this work, we investigate the impact of the non-resonant continuum on 
single-nucleon overlap integrals and spectroscopic factors in weakly bound 
and unbound nuclear states. Our paper is organized as follows. The model 
used is the SM in the complex momentum-plane, the so-called Gamow Shell 
Model (GSM) [7,8,9], which is briefly presented in Sec. 2. Section 3 describes 
the concrete realization of the GSM employed in this work, i.e., interactions 
used and the GSM model space, as well as the approximate schemes that 
have been introduced to illuminate speciflc physics points. The discussion of 
spectroscopic factors and overlap integrals is contained in Sec. 4. We show that 
the realistic radial overlap functions, which are needed for the description of 
transfer reactions, can be generated by single-particle wave functions of the 
appropriately chosen complex average potential reproducing the complex Q- 
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value of the studied reaction or decay. Furthermore, we demonstrate that the 
scattering continuum is of crucial importance for the spectroscopic factors in 
the vicinity of the particle emission threshold. This problem has been discussed 
recently in our work [10] for the i = 1 neutron waves. Here, we extend this 
discussion to £ = 2 and also compare the spectroscopic factors in GSM with 
those calculated in the pole approximation to see the influence of non-resonant 
continuum. As an illustrative example in these studies, we choose the case of 
model two-neutron systems outside the inert core: ^He and ^^O. Our aim is 
not to fit the actual experimental data, but rather to illustrate generic effects. 
Finally, the conclusions of this work are summarized in Sec. 5. 



2 Gamow Shell Model 



As the Gamow Shell Model has been described in detail in a number of publica- 
tions, only rudimentary information is provided here. In the roots of GSM lies 
the Berggren one-body completeness relation [11,12] that provides the math- 
ematical foundation for unifying bound and unbound states. The Berggren 
ensemble consists of bound, resonant, and scattering s.p. wave functions, gen- 
erated by a finite-depth potential V{r), either real or complex. The wave 
functions are regular solutions of the s.p. Schrodinger equation, 

that obey outgoing or scattering boundary conditions: 

UB{r)r^+oo = C+ H+{kr) + C_ Hi{kr), (2) 




where k = \/2mE/h is the complex wave number. For the resonant states, 
C_=0 (outgoing boundary condition). In this paper, as we consider only va- 
lence neutrons, Hf" is a Hankel function. Normalization constants C_|_ and C_ 
are determined from the condition that the radial wave functions u{r) are 
normalized to unity (for resonant states) or to the Dirac delta (for scattering 
states). 

For a given partial wave the scattering states are distributed along the 

contour in the complex momentum plane. The set of all bound and 
resonant states enclosed between and the real /c-axis, and scattering states 
is complete [11]: 

Uts){uB\ = l. (3) 
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The Gamow (and Berggren) states are vectors in the rigged Hilbert space (or 
Gelfand triplet) [13,14,15,16]. Consequently, in the above, the tilde symbol 
signifies that the complex conjugation arising in the dual space affects only 
the angular part and leaves the radial wave function unchanged [11,12]. The 
Berggren states are normalized using the squared radial wave function and 
not the modulus of the squared radial wave function. In the standard SM 
treatment, the latter normalization is used. 

In numerical applications, the integral over scattering states along L]_ has to 
be discretized and the selected scattering states have to be renormalized [7]. 
This leads to a discrete completeness relation: 

N 

^|uB)(n^|~l ; \ub) = ^f^W^), (4) 

where {ks^'^B) is the set of discretized complex wave numbers and associated 
weights provided by a Gauss-Legendre quadrature. As discussed in Ref. [17], 
the Gauss-Legendre integration formula offers an excellent precision at a mod- 
est number (~30) of discretization points. 

The many-body Berggren basis is that of Slater determinants built from s.p. 
states of Eqs. (2-4): 

\SD.>) = \ui^ ■ --Ui^), (5) 

where the index i labels the many-body basis, and Ui- is the j-th s.p. state 
occupied in \SDi). The many-body completeness relation is a consequence of 
Eq. (5); it is obtained by forming all possible many-body Slater determinants: 

Y,\SD,){SD,\c^l. (6) 

i 

The equality in Eq. (6) is not exact due to the discretization of the one-body 
completeness relation. 

In the Berggren representation, the SM Hamiltonian matrix H becomes com- 
plex symmetric. Its complex-energy eigenvectors can be obtained by using the 
complex Lanczos method [18]. The fundamental difference between GSM and 
the real-energy SM is that the many-body resonant states of the GSM are 
embedded in the background of scattering eigenstates, so that one needs a cri- 
terion to isolate them. The overlap method introduced in Ref. [7] has proven 
to be very efficient in this respect. To this end, one diagonalizes H in the basis 
spanned on the s.p. resonant states only (the so-called pole approximation); 
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this yields a zeroth-order approximation |\l/o) for the wave function. The vec- 
tor |\l/o) has a correct outgoing behavior; it is used as a pivot to generate a 
Lanczos subspace of the full GSM space. The requested resonant eigenstate of 
H is the one which maximizes the overlap |(\l/o|^')|. The overlap method has 
also been employed in the Density Matrix Renormalization Group (DMRG) 
technique [19] recently generalized to treat the non-hermitian GSM case [20]. 

The definition of observables in GSM follows directly from the mathematical 
setting of quantum mechanics in the rigged Hilbert space [13,14,15,16] rather 
than the usual Hilbert space. A modified definition of the dual space, embodied 
by the tilde symbol above the bra states in Eqs. (3-6), implies that observ- 
ables in many-body resonances become complex. In this case, the real part 
of a matrix element corresponds to the expectation value, and the imaginary 
part can be interpreted as the uncertainty in the determination of this expec- 
tation value due to the possibility of decay of the state during the measuring 
process [11,16,21]. That is a standard price for replacing the time-dependent 
description of an unstable system by the quasi- stationary description. 



3 Gamow Shell Model Implementation 

This section contains details pertaining to GSM calculations carried out in this 
work. It describes the Hamiltonian used, the choice of the active Hilbert space, 
and the approximations employed. All remaining details of the formalism can 
be found in Refs. [7,10]. 



3. 1 GSM Hamiltonian 

The doubly magic nuclei ^He and ^^O are assumed to be the inert cores in 
our helium and oxygen calculations, respectively. The GSM Hamiltonian con- 
sists of a one-body Woods-Saxon (WS) potential, representing a core-valence 
interaction, and a two-body residual interaction: 



The WS potential contains the central term and the spin-orbit term: 




(7) 



V{r) = -Vws- f{r)-AVso (l-s) 



1 df\r) 
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r dr 
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where Vws and Vso are the strength constants and 



fir) 



1 + exp 



d 



(9) 



is the WS form-factor characterized by radius Ro and diffuseness d. In the 
hehum calculations, we employed the "^He" WS parameter set [7] which re- 
produces the experimental energies and widths of known s.p. resonances 3/2j~ 
and l/2i in ^He {Op3/2 and Opi/2 resonant states in our model). In the oxygen 
case, we used the "^'''O" WS parameter set [7] that reproduces experimental 
energies and widths of the 5/2i and l/2i bound states and the 3/2f resonance 
in -^'^O (i.e., 0^5/2, lsi/2, and 0^3/2 resonant states). 

We used a different residual interaction for helium and oxygen isotopes. This 
choice has been motivated by our previous studies of oxygen, helium, and 
lithium chains [7,8]. For the heliums, we took the finite-range surface Gaussian 
interaction (SGI) [8] 



m 



2Ro) 



(10) 



with the range /i=lfm and the coupling constants depending on the total 
angular momentum J of the neutron pair: 

—392 MeV fm^. These constants are fitted to reproduce the ground-state (g.s.) 
binding energies of ^He and ^He in GSM. 

For the oxygens, we employed a surface delta interaction (SDI) [8]: 

V}!I = E H ■ ^ (rr - rS) ■ 5(|rl| - Ro), (11) 

i<j 



where the SDI coupling constant Vo=-700 MeV fm'^ was fitted to in order to 
reproduce the two-neutron separation energy of ^^O. 



3.2 GSM Configuration Space 



In the helium variant, the valence space consists of the P3/2 and pi/2 neu- 
tron partial waves. The ps/2 wave functions include a 0^3/2 resonant state 
and P3/2 non- resonant scattering states along a complex contour enclosing the 
0p3/2 resonance in the complex /c-plane. For a pi/2 part, we take non-resonant 
scattering states along the leal-k axis (the broad Opi/2 resonant state plays a 
negligible role in the g.s. wave function of ^He). For both contours, the maxi- 
mal momentum value is /c^ax =3.27 fm~^ The contours have been discretized 
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with up to 60 points, and the attained precision on energies and widths is 
better than 0.1 keV. 

In the oxygen case, the valence GSM space for neutrons consists of the 0^5/2, 
lsi/2, and 0(^3/2 Gamow states and the non-resonant complex 6/3/2 continuum. 
The 0d5/2 and lsi/2 shells are bound while the 0^3/2 state is a narrow res- 
onance. The maximum value for k on the contour L'^^^ is fcmax =1-5 fm~^. 
The contour has been discretized with 45 points and the resulting precision 
on energies and widths is ~ 0.1 keV. The real-energy ^5/2 and Si/2 contours 
give rise to a renormalization of the residual interaction coupling constant; 
hence, they have been neglected [7]. 

3.3 Simplified SM schemes: GSM-p and HO-SM 

In order to illuminate certain aspects of GSM results (e.g., non-resonant con- 
tinuum coupling or configuration mixing), we have introduced two simplified 
SM schemes: the GSM in the pole approximation (GSM-p) and the harmonic 
oscillator shell model (HO-SM). 

In GSM-p, the scattering components of the Berggren basis are disregarded, 
and only the resonant states are present in the basis. In particular, they will 
be the OP3/2 and Opi/2 neutron states for heliums, and Oo?5/2, lsi/2, and 0^3/2 
neutron states for oxygens. In this case, the one-body completeness relation (3) 
is obviously violated. Still, the comparison between full GSM results and GSM- 
p is instructive as it illustrates the influence of the non-resonant continuum 
subspace on the GSM results. 

For the HO-SM calculation, the radial wave functions of the GSM-p basis are 
replaced by those of the spherical harmonic oscillator (HO) with the frequency 
hu = 41^-^/3 MeV. The real energies of the resonant states define the one- 
body part of the HO-SM Hamiltonian. The HO-SM scheme is intended to 
illustrate the "standard" CQS SM calculations in which only bound valence 
shells are considered in the s.p. basis. 



4 One-nucleon overlap integrals and spectroscopic factors in GSM 

4-1 Definitions 

Single-nucleon overlap integrals and the associated spectroscopic factors are 
basic ingredients of the theory of direct reactions (single-nucleon transfer, nu- 
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cleon knockout, elastic break-up) [22,23,24,25]. Experimentally, spectroscopic 
factors can be deduced from measured cross sections; they are useful mea- 
sures of the configuration mixing in the many-body wave function. The asso- 
ciated reaction-theoretical analysis often reveals model- and probe- dependence 
[26,27,28] raising concerns about the accuracy of experimental determination 
of spectroscopic factors. In our study we discuss the uncertainty in deter- 
mining spectroscopic factors due to the two assumptions commonly used in 
the standard SM studies, namely (i) that a nucleon is transferred to/from a 
specific s.p. orbit (corresponding to an observed s.p. state), and (ii) that the 
transfer to/from the continuum of non-resonant scattering states can be dis- 
regarded. This discussion complements our recent work [10] whose focus was 
on threshold anomalies and channel coupling. 

The one-nucleon radial overlap integral is discussed in the usual way [23,24,25]: 



u,,{r) = {^'/{n,---,rA = f^\ |vl>;;i-i^(fi,---,rVi))® "1. (12) 



.Ja 



Ja, 



where l^^-^) and I^X."/) are wave functions of nuclei A and A — 1, respectively, 
and angular and spin coordinates (represented by are integrated out so 

that U£j depends only on r = \ f\. The spectroscopic factor is given by the real 
part of the norm S"^ of the overlap integral. Using a decomposition of the s.p. 

channel in the complete Berggren basis, one obtains GSM expressions 
for the overlap integral and S"^: 

=^{<'\\4m\<"-i) ^^^3\ub) (13) 
S'^j ul{r)dr =^{^/\\a%{BmA-if (14) 

where al-{B) is a creation operator associated with a s.p. basis state \uq). 
We wish to emphasize that since Eqs. (13,14) involve summation over all 
discrete Gamow states and integration over all scattering states along the 
contour L^, the final result is independent of the s.p. basis assumed (see 
also Ref. [29] where a model-independent spectroscopic factor experimental 
extraction procedure is proposed). This is in contrast to the standard SM 
treatment of spectroscopic factors where model-dependence enters through 
the specific choice of a s.p. state [31], with Eq. (14) reducing to a single 

matrix element \{'^^a\\ o-tij 1 1 ^ a-Y ) P • 

Spectroscopic factors are often extracted from measured transfer reaction cross 
sections. (In the context of the following discussion, is worth noting that while 
the reaction cross sections are measurable quantities, the spectroscopic factors 
are purely theoretical concepts. First, they are deduced from experiment in 
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a model-dependent way. Second, occupation numbers deduced from spectro- 
scopic factors are not observables per se [30].) In the CQS SM framework, the 
transfer total cross section is given by [32,24,25]: 

= H Snijm (15) 
nijm 

where Snijm is the spectroscopic factor associated with the \n£jm) s.p. state, 
and cT^p^"* a s.p. cross section containing all the kinetic dependence related to 
the transfer of the nucleon of the projectile to the shell model state. While 
Eq. (15) is, in principle, exact for bound s.p. states, as one sums over all 
quantum numbers nijm, the commonly used approximation is to take only 
one radial quantum number n for each ijm in the sum of Eq. (15), i.e., one 
assumes that the nucleon is transferred to/from an orbit with a specific radial 
quantum number. Hence, the resulting value of spectroscopic factor becomes 
spuriously basis-dependent. We shall show in the following that this can lead 
to sizeable errors if the states l'^'^^) or I^I/^'l"/) are loosely bound. 

Obviously, the factorization (15) does not apply to many-body states which 
cannot be described in the CQS formalism. Firstly, in the OQS formalism the 
sum in Eq. (15) is not limited to discrete states \nijm) only. The appearance 
of the cusp at the channel threshold both in the reaction cross-section and in 
the real part of the GSM spectroscopic factor suggests an expression: 

a = J2ReiSnej)a%,, (16) 

as the relation between the observable reaction cross-section and the real part 
of the spectroscopic amplitude S = Snij (cf (14)) calculated in GSM. Secondly, 
as stated earlier, observables involving unbound states are complex in GSM. 
It was conjectured by Berggren [11] and later by Gyarmati et al. [33] that 
the real part of a matrix element corresponds to the expectation value. Based 
on scattering theory, Berggren showed [34] that the real part of the complex 
cross section for populating a resonance is equal to the energy integral of the 
in-elastic continuum cross section across the resonance peak. The imaginary 
part of the cross section can be identified with the strength of the resonance- 
background interference. This conclusion was later generalized to hold for 
any observable operator involving resonant states. In this work, we follow 
Berggren's interpretation and associate an imaginary part of a matrix element 
(e.g., spectroscopic factors or radial overlap function) with the uncertainty in 
the determination of the real part of this matrix element. 

For weakly bound or unbound nuclei, the transfer reaction cross section is 
changed by couplings to the non-resonant continuum which modify both the ef- 
fective interaction among valence particles and the configuration mixing. This 
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change in the cross section has a regular component having a smooth depen- 
dence on the distance from the particle-emission threshold, and an anomalous 
component which results in the Wigner-cusp phenomenon [35]. Other quanti- 
ties, such as the strength function [36,37], the continuum-coupling correction 
to the CQS eigenvalue [38], or the overlap integral [10] are modified in a sim- 
ilar way. These two components can be seen by comparing full GSM results 
with equivalent SM results (here: the HO-SM approximation). 

Generic features of wave functions in the vicinity of the reaction threshold 
(the particle-emission threshold) are responsible for the anomalous component. 
Wigner was first to notice it as the characteristic behavior of scattering and 
reaction cross-sections [35] , often referred to as the Wigner threshold law or the 
Wigner-cusp phenomenon. A quantitative explanation of this behavior of cross 
sections was given later [39,40,41,42,43,44]. Wigner threshold law tells that in 
the vicinity of particle-emission thresholds the reaction cross section has a 
universal shape depending on the charge and the orbital angular momentum 
of the scattered particle. In particular, the reaction cross section for neutron 
scattering near the neutron emission threshold shows a characteristic cusp 
dependence on the orbital momentum of a scattered neutron. This salient 
dependent anomalous behavior (cusp) is absent for charged scattered particles. 

Contrary to the universality of the cusp, its magnitude is not a generic feature 
of OQSs and cannot be found from general principles. Also the ratio between 
regular and anomalous components in the near-threshold behavior of reac- 
tion cross-sections, strength functions, radial overlap integrals, spectroscopic 
factors, or continuum-coupling correction to CQS eigenvalues is the specific 
property of each considered OQS and its Hamiltonian. 

The modification of the cross section in one channel manifests itself due to 
the unitarity of the scattering matrix in other open channels. This coupling 
effect has been studied theoretically [40,41,42,43,45,36,10] and experimentally 
[46,47,48,49] in connection with the Wigner threshold anomaly. The channel- 
channel continuum coupling influences both the anomalous and regular cross- 
section components. Similar multi-channel threshold effects are expected to 
show up in various observables which depend on occupation probabilities of 
s.p. shells and the configuration mixing involving continuum states. 

It was shown in Ref. [10] that once a given nuclear Hamiltonian is chosen, 
quantities such as the spectroscopic factors, defined in GSM through the norm 
of the overlap integral, and their behavior close to particle-emission thresholds 
are defined uniquely in terms of the exact many-body GSM solutions. Below, 
we shall demonstrate that the threshold variations due to the continuum cou- 
pling can be of a comparable size to those generated by the long-range corre- 
lations and the coupling to high-momentum states reached by short-range and 
tensor components of the nucleon-nucleon interaction [50]. Hence, the states 
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e (MeV) 

^3/2 

Fig. 1. The real part of the overlap (^He(g.s.) | [^He(g.s.) ^Vz/'i^^^)'^ ^ ^ function of 
the energy of the 0^3/2 resonant state. ^He is always bound. The solid line (GSM) 
shows the full GSM result. The dotted line (HO-SM) corresponds to the SM cal- 
culation in the oscillator basis of 0p3/2,0pi/2- The dashed line (GSM-p) shows the 
GSM result in the pole approximation. The imaginary part of S*^ is given in the 
inset. In order to better illustrate the effect of configuration mixing, S"^ has been 
normalized to 1 in the limit of vanishing residual interaction. 

of open and closed quantum systems may have not only different asymptotic 
behavior but also different shell model structure. In that sense, loosely bound 
or unbound states form a different class of correlated quantum many-body 
systems than those found in well-bound nuclear states. 

4.2 The ^He+n^^He case 

We first consider the g.s. spectroscopic factor of ^He corresponding to the 
channel [^He(g.s.) ^^3/2]°^- In order to investigate the effect of the continuum 
coupling, the depth of the WS potential is varied so that the energy Cp^^^ of 
the 0p3/2 s.p. state (the lowest ^3/2 pole of the S'-matrix) changes its character 
from bound to unbound. In our model space, the energy e^g^^ both the g.s. 
energy of ^He and (negative of) the one-neutron (In) separation energy Sin 
of ^He. 

The calculated spectroscopic factor pHe(g.s.)®p3/2]°^ in ^He(g.s.) is shown in 
Fig. 1 as a function of the real energy Cp^^^ of the 0^3/2 pole. In the range of 
^P3/2 values shown in Fig. 1, ^He is bound (S'i„(^He)>0). The spectroscopic 
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factor strongly depends on the position of the pole: for ep._^^^<0 (S'i„(^He)>0) 
it decreases with Cp^^^ while it increases for ep^^^>0 (S'i„(^He)<0). At the In- 
emission threshold in ^He (ep3^2=0), spectroscopic factor exhibits a cusp. At 
this point, the derivative of spectroscopic factor becomes discontinuous and 
the coupling matrix element between the resonant 0^3/2 state and the non- 
resonant continuum reaches its maximum [38]. 

As discussed in Ref. [10], the dependence of spectroscopic factor on Cp^^^ 
around eop3/2=0 (where the spectroscopic factor exhibits a cusp due to a 
coupling with the ^He+n+n channel) follows the threshold behavior of the 
reaction cross section [35]. Specifically, the anomalous component of spectro- 
scopic factor below the In threshold in ^He behaves as (— e^j)^"^/^. Above the 
threshold, the spectroscopic factor factor is complex; the real part behaves 
as (e^j)^"*"^/^ while the imaginary part, associated with the decaying nature of 
^He, varies as {eijY"^^"^. 

To assess the role of the continuum, both resonant and non-resonant, the GSM 
results are compared in Fig. 1 with the HO-SM and GSM-p calculations. In 
contrast to GSM, spectroscopic factors in HO-SM and GSM-p vary little in 
the energy range considered and no threshold effect is seen. It is interesting 
to note that in the limit of an appreciable binding, the 0^3/2 wave function 
is fairly well localized, the importance of the continuum coupling is dimin- 
ished, and the HO-SM result approaches the GSM limit. For Cp^^^ > 0, GSM 
results are in strong disagreement with both HO-SM and GSM-p. A ~25% 
difference between spectroscopic factors in HO-SM and GSM-p is striking. 
This difference reflects a strong dependence of the two-body matrix elements 
on the actual radial wave functions. In the GSM-p variant, the g.s. of ^He 
is described by an almost pure [0p3/2 ® 0^3/2]°^ configuration and the result- 
ing spectroscopic factor is close to one in the whole energy region considered. 
A difference between GSM and GSM-p results illustrates the impact of the 
non-resonant continuum. 

Figures 2-4 display one-neutron radial overlap integrals calculated in GSM 
for bound ^He and three energies of the 0^3/2 pole corresponding to differ- 
ent physical situations considered in Fig. 1. In Fig. 2, both ^He and ^He are 
bound with one-neutron separation energies 0.5 MeV and 1 MeV, respectively. 
Figure 3 illustrates the case of unbound ^He (-E[^He]=0.25 -i0.056 MeV) and 
bound ^He (-E[^He]=-0.5 MeV). Here, the generalized one-neutron separation 
energy. Sin, defined as a difference between g.s. binding energies of neighboring 
nuclei, becomes complex: 

~Sin{N) = E{N - 1) - E{N) = SUN) - [r{N - 1) - T{N)] . (17) 
The real part of is the usual separation energy Sin while the imaginary 
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Fig. 2. Single- neutron 
(6He(g.s.)|[5He(g.s.) ® ^3/2]°' 



radial overlap integral n(r) for the 
^ channel calculated in GSM for one-neutron 
bound ^He and ^He. The radial wave function of the 0^3/2 resonant state of the 
real WS potential adjusted to reproduce the GSM value of Sin in ^He practically 
coincides with u{r). The overlap integral has been normalized to unity to allow 
comparison with the WS wave function. 

part appears when either parent or daughter nucleus is unbound. Finally, Fig. 
4 shows u{r) for the case of bound ^He and ^He having the same energy, i.e.. 

The overlap integral is peaked around the surface (~2 fm) and its asymptotic 
behavior is determined by one- neutron separation energy [51,52]: 

u{r) e-^^ (18) 



where the decay constant is 



K = J2mSin/h. 



(19) 



In the case considered, u{r) is expected to decay exponentially with k deter- 
mined by the one-neutron separation energy of ^He. Figures 2-4 show, by the 
dotted line, the radial wave function of the 0p3/2 resonant state of the WS 
potential with a depth adjusted to reproduce the GSM value of Sin in ^He. 
Since ^He is one-neutron bound, the corresponding 0^3/2 WS state is bound; 
hence, its wave function is real. It is seen that the agreement between the WS 
wave function and the real part of u{r) is excellent. 

The situation shown in Fig. 3 is particularly interesting as the overlap integral 
is complex due to an unbound nature of ^He. Here, a real WS potential cannot 
provide any information about the imaginary part of u{r). One possible way of 
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Fig. 3. Solid line: the single-neutron radial overlap integral u{r) for the 
(^He(g.s.)l [^He(g.s.) ^ P3/2]^ )^ channel calculated in GSM for unbound ^He and 
bound ^He. Dotted line: the radial wave function of the 0^3/2 resonant state of 
the real WS potential whose strength was adjusted to reproduce the GSM value 
of Sin in ^He. Dashed line: the radial wave function of the 0^3/2 resonant state of 
the complex WS potential whose complex strength was adjusted to reproduce the 
GSM value of the generalized separation energy in ^He. All wave functions are 
normalized to unity to allow comparison. The real (imaginary) parts of the wave 
function are shown in the top (bottom) panel. 
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Fig. 4. Similar to Fig. 2 except for the threshold situation of 5'i„pHe]=0. 

overcoming this difficulty is to introduce a complex WS potential characterized 
by a complex strength: 



Vws = v^s + 



WS- 



(20) 



Following the usual SM strategy applied to real u{r) cases, Vws can be ad- 
justed to reproduce the generalized separation energy (17). This guarantees 
that the asymptotic behavior of the overlap integral, given by Eqs. (18-19) with 
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E[^He] - E[^He] (MeV) 

Fig. 5. Similar to Fig. 1 except as a function of (negative) Sin of ^He. The depth of 
the WS potential of ^He has been adjusted to bind ^He (see text). The separation 
energy Sin of ^He has been varied by changing the coupling constant Vq'^~^^ of the 
SGI two-body interaction. 

Sin-^Sin, is correct. (The extension of the discussion of asymptotic properties 
of u{r) in Ref. [51] to the complex-energy case is straightforward; it follows 
from the analytic properties of the vertex form factor.) It is seen in Fig. 3 
(dashed line) that both real and imaginary parts of the radial overlap integral 
are well reproduced by the 0^3/2 resonant wave function of a WS potential 
that reproduces the GSM value of ^in. This 0^3/2 complex-energy state cor- 
responds to a 'decaying bound state' which can be found only in the complex 
s.p. potential (see, e.g., [53,54]); asymptotically, such a state exhibits an ex- 
ponentially damped oscillation. The effect of the non-resonant continuum is 
seen in a slightly better localization of the GSM overlap. 

Figure 4 illustrates the threshold limit. Here, 5'i„=0 and u{r) shows asymptotic 
behavior that is not exponential. Also in this case the GSM radial overlap 
integral is well reproduced by the radial 0^3/2 WS wave function corresponding 
to a halo state located at zero energy. 

The Wigner estimate for the near-threshold behavior of cross section [35] is 
independent of the reaction mechanism. While the Wigner limit has been 
reached in Fig. 1 by changing the pole of the one-body S-matrix, it is in- 
teresting to see whether a similar pattern can be generated by many-body 
correlations. Figure 5 illustrates a direct case of the Wigner cusp. Here, we 
have fixed the WS potential so that the 0p3/2 and Opi/2 shells are both bound 
with respective energies of —5 MeV and —0.255 MeV, and we have varied the 
SGI coupling constant Vq'^~^^ so that 5*1,1 [^He] changes sign. The behavior of 
spectroscopic factor around the In threshold is similar to that of Fig. 1, with 
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the GSM spectroscopic factor exhibiting a non-analytic behavior at S'i„=0. 
This result constitutes an excellent test of the GSM formalism: the Wigner 
limit is reached precisely at a threshold obtained from many-body calcula- 
tions. In the whole range of considered separation energies the HO-SM and 
GSM-p results vary smoothly and they differ from the GSM results. In this ex- 
ample, the anomalous component of spectroscopic factor is much smaller than 
the regular component showing the essential role played by the non-resonant 
continuum component in the ground state of ^He. 




r(fm) 

Fig. 6. Similar to Fig. 3 except for bound ^He and unbound ^He {Sin of ^He is 
negative.) The dashed line depicts the radial wave function of the 0^3/2 resonant 
state of the complex WS potential whose depth is adjusted to reproduce the GSM 
value of the generalized separation energy Sin in ^He. 

Figures 6 and 7 display u{r) calculated for the two cases shown in Fig. 5 repre- 
senting unbound ^He, and ^He being either bound (Fig. 6) or unbound (Fig. 7). 
In Fig. 6, the generalized separation energy is complex, S'i„=-0.2-|-i0.0915 MeV. 
The corresponding Gamow pole is a bound decaying state. Figure 7 shows the 
radial overlap integral for unbound ^He and ^He having the same real ener- 
gies, i.e., Sin=0, but different widths. Consequently, the generalized separation 
energy Sin has a non-zero imaginary part in this case. The corresponding res- 
onant state lies on the diagonal of the fourth quarter of the complex-/c plane. 

The radial wave functions of the 0^3/2 resonant state of the real WS potential 
adjusted to the GSM separation energy of ®He are shown by the dotted lines in 
Figs. 6 and 7. The agreement between the WS wave function and u{r) is poor. 
In particular, their real parts have a different asymptotic behavior. In Fig. 
7, the imaginary part of the 0^3/2 WS wave function is zero by construction, 
and it does not reproduce the salient behavior of the imaginary part of u{r). 
On the other hand, the In radial overlap integral is very well reproduced by 
the 0p3/2 resonant state wave function of the complex WS potential (dashed 
curves in Figs. 6, 7). 
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Fig. 7. Similar to Figs. 3 and 6 except for unbound ^He and ^He. The In separation 
energy Sin of ^He has been adjusted to zero. The generahzed In separation energy 
Sin has a non- vanishing imaginary part. 

In general, we find that shapes of In overlap integrals calculated in GSM 
for any complex value of Sin are well reproduced by radial wave functions of 
a complex WS potential whose depth is adjusted to reproduce the complex 
separation energy. In particular, for unbound parent and/or daughter nuclei, 
the asymptotic behavior of u{r) is given by the generalized complex separation 
energy Sin and not by the real separation energy Sin- 



4.3 The^'^O+n^^^O case 



To investigate the dependence of spectroscopic factors on the orbital angular 
momentum £ of the transferred nucleon, we consider the ^3/2 partial wave in 
^®0. We show in Fig. 8 the spectroscopic factor for the excited Ojj' state of 
^^O in the channel [^'''0(3/2^) ® ^3/2]''^. The behavior of spectroscopic factor 
shown in Fig. 8 is similar to that of Fig. 1, except the variations are much 
weaker (cf. the dramatically expanded scale) and the threshold behavior is 
different. Namely, the spectroscopic factor is continuous and smooth; it is 
its derivative that exhibits a cusp around the 0^3/2 threshold. Again, this 
is consistent with the general expectation that, for £=2, Re{S'^) below the 
In threshold of ^'^O behaves as (— e^j)^/^ while above the threshold Re{S'^) 
(/m(S'^)) should behave as (e^j)^/^ {{^ij)^^'^)- The associated In ^3/2 radial 
overlap integrals (not displayed) are extremely close to the s.p. resonant 0^3/2 
wave functions. 

To investigate the behavior of spectroscopic factor in ^^O around S'i„=0, the 
depth of the WS potential has been decreased to lower the position of the 
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e^,.. (MeV) 



'3/2 



Fig. 8. Similar to Fig. 1 except for the overlap in the excited 0^ state of ^^O: 
(^^O(0;[)|[^'^O(3/2^") (g) 4/2]°^)^. The first derivative of the spectroscopic factor in 
the neighborhood of the ed^^^=0 threshold is shown in the inset (a) while inset (b) 
displays the imaginary part of 5^. 




0.05 0.05 

A (MeV) 



E[i^O*] - £[^"^0*] (MeV) 



Fig. 9. Similar to Fig. 5 except for the overlap in the excited O3 state of O: 
(^^^0(03^)1 [^^"^0(3/2^) (g) d^/2f^)^- The energy derivative of the spectroscopic factor 
around the In threshold is shown in the inset. 

Od3/2 s.p. pole down to -5 MeV. At the same time, the SDI coupling strength 
Vq has been modified to allow Sin of the O3 state of ^^O to go through zero. 
The results displayed in Fig. 9 are consistent with the Wigner limit for i=2. 
Namely, spectroscopic factor and its first derivative are continuous around 
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5'in=0, while the second derivative exhibits discontinuity (see the inset in 
Fig. 8). The associated In 0^3/2 radial overlap integrals (not displayed) are 
extremely close to the s.p. resonant 0^3/2 wave function of a real (or complex) 
WS potential. 



5 Conclusions 

By explicit many-body GSM calculations that fully account for a coupling to 
the scattering space, we demonstrate the presence of a near-threshold, non- 
perturbative rearrangement in the wave function that has an appreciable low 
angular momentum s.p. component. The threshold behavior of spectroscopic 
factors (Wigner cusp) can only be reproduced through a complete inclusion of 
scattering states, including the non-resonant space, in the GSM basis. Having 
a complete basis which allows to describe bound, weakly-bound and unbound 
states on the same footing, is the only way to guarantee the unitarity, which 
lies at the basis of the Wigner threshold effect in cross sections and other 
observables, in particular in the multichannel case. This fundamental require- 
ment is not respected in any CQS formulation of the many-body theory. As 
shown in several examples discussed in our work, restoration of unitarity can 
strongly affect both values and behavior of spectroscopic factors. 

The detailed analysis of near-threshold behavior of spectroscopic factors asso- 
ciated with £=1 (heliums) and i=2 (oxygens) partial waves shows the presence 
of the Wigner limit in the many-body GSM solution. Namely, the fluctuat- 
ing component of the spectroscopic factor behaves as {—eejY~^^'^ below the 
In threshold. Above the threshold, the spectroscopic factor is complex; the 
real part behaves as (e^jY^^^'^ while the imaginary part, associated with the 
decaying nature of a resonance, varies as (e^j)^"^/^. 

If either parent or daughter nucleus is unbound, the corresponding one-nucleon 
overlap integral is complex. It can be very well approximated by a resonant 
state of a one-body potential which reproduces the complex generalized one- 
nucleon separation energy Sin (17). This is a straightforward generalization of 
an approximate treatment of one-nucleon overlaps often used in SM studies 
(see, e.g., Ref. [28]) where a real average potential is employed with a depth 
adjusted to reproduce Sin- Therefore, we conclude that the realistic radial 
overlap functions, which constitute the basic theoretical input needed for a 
description of transfer reactions, can be conveniently generated by a radial 
s.p. wave function of the (complex) one-body potential reproducing the gener- 
alized separation energy. On the other hand, the normalization of the overlap 
integral, i.e., the spectroscopic factor, cannot be obtained in a simple way; 
here the full microscopic treatment is necessary (see discussion in Ref. [10]). 
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Optimally, the complex average potential used to generate one-nucleon overlap 
functions should be generated self-consistently using the Gamow-Hartree-Fock 
procedure [8] . The imaginary part of this potential should not be confused with 
the absorption potential used in the context of optical model studies. As we 
illustrated in this paper, depending on Sin, the corresponding resonant states 
do not always obey the usual Berggren classification for hermitian Hamiltoni- 
ans. For instance, one has to consider decaying bound states (resonant states 
lying in the first quarter of the complex-A; plane) or threshold states that have 
zero energy but nonzero width. 

While Berggren never considered complex potentials in the derivation of com- 
pleteness relation for Gamow states [11,12], its demonstration in the case of 
localized potentials, as those of the WS type used in this paper, is straightfor- 
ward [55]. For complex potentials, multiple S*- matrix poles can appear in the 
resonant-state expansions. This complicates the completeness relation as some 
additional states, which are not eigenf unctions of the one-body Hamiltonian, 
have to be added to reach completeness [55]. Moreover, a possible appearance 
of the spectral singularities, i.e., the S'-matrix poles lying on the real momen- 
tum axis, not considered in Ref. [55], can complicate matters [56]. However, 
multiple singularities cannot happen if the imaginary part of the potential is 
small enough, and - in this case - the Berggren completeness relation for com- 
plex potentials is obtained in the same manner as for real potentials. Such a 
relation could be useful when developing the generalized Gamow-Hartree-Fock 
method that deals with unusual Gamow states. 
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